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Hydrogen wall and heliosheath Ly-a absorption toward nearby 
stars: possible constraints on the heliospheric interface plasma 
flow 
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Abstract. In this paper, we study heliospheric Ly-a absorption toward nearby stars in dif- 
ferent lines of sight. We use the Baranov-Malama model of the solar wind interaction with 
a two-component (charged component and H atoms) interstellar medium. Interstellar atoms are 
described kinetically in the model. The code allows us to separate the heliospheric absorp- 
tion into two components, produced by H atoms originating in the hydrogen wall and heliosheath 
regions, respectively. We study the sensitivity of the heliospheric absorption to the assumed 
interstellar proton and H atom number densities. These theoretical results are compared with 
interstellar absorption toward six nearby stars observed by the Hubble Space Telescope. 



1. Introduction 

The solar system is traveling in the suiTounding Local Interstel- 
lar Cloud (LIC). In the 1960s, it was realized [e.g., Patterson et 
ah, 1963; Fahr, 1968; Blum andFahr, 1970] that interstellar atoms 
penetrate deep into the heliosphere and, therefore, can be observed. 
The interstellar atoms of hydrogen have been detected by measure- 
ments of the solar backscattered Ly-a irradiance [Bertaux and Bla- 
mont, 1971; Thomas and Krassa, 1971]. Later, interstellar atoms 
of helium were also measured directly [Witte et al, 1996] and in- 
directly as backscattered solar irradiance [Weller and Meier, 1981; 
Dalaudier et al, 1984]. At present, there is no doubt that inside the 
heliosphere, properties of He atoms such as temperature and veloc- 
ity are different from those of H atoms. In particular, the H atoms 
are decelerated and heated compared with atoms of interstellar he- 
lium inside the heliosphere [Lallement et al., 1993; Costa et al, 
1999; Lallement, 1999]. The reason for this is a stronger coupling 
of H atoms with plasma protons in the heliospheric plasma interface 
through charge exchange. Many works developed the concept of 
the heliospheric plasma interface over more than four decades after 
pioneering papers by Parker [1961] and Baranov et al. [1970]. 

The heliospheric interface is formed by the interaction of the 
solar wind with the charged component of the interstellar medium 
(see Figure 1). The heliospheric interface is a complex structure 
having a multi-component nature. The solar wind and interstellar 
plasmas, interplanetary and interstellar magnetic fields, interstel- 
lar atoms. Galactic and anomalous cosmic rays (GCRs and ACRs), 
and pickup ions all play prominent roles. Interstellar hydrogen 
atoms interact with the heliospheric interface plasma by charge ex- 
change. This interaction significantly influences both the structure 
of the heliospheric plasma interface and the flow of the interstellar H 
atoms. In the heliospheric interface, atoms newly created by charge 
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exchange have the properties of local protons. Since the plasma 
properties are different in the four regions of the heliospheric in- 
terface shown in Figure 1, the H atoms in the heliosphere can be 
separated into four populations, each having significantly different 
properties. For example, population 3 consists of the atoms created 
by charge exchange with relatively hot protons in the region of dis- 
turbed interstellar plasma around the heliopause (region 3 in Figure 
1). It was realized theoretically by Baranov et al. [1991] that the 
atoms of population 3 form a so-called "hydrogen wall" around the 
heliopause. The hydrogen wall is a significant enhancement of the 
density of interstellar atoms of hydrogen around the heliopause com- 
pared with the number density of the undisturbed local interstellar 
medium. The self-consistent models of the heliospheric interface 
predict that the secondary interstellar atoms (or atoms of population 
3) are decelerated and heated compared with the original interstellar 
H atoms. 

The Ly-a transition of atomic H is the strongest absorption line 
in stellar spectra. Thus, the heated and decelerated atomic hydro- 
gen within the heliosphere produces a substantial amount of Ly-a 
absorption. This absorption, which must be present in all stellar 
spectra since all lines of sight go through the heliosphere, has been 
unrecognized until recently, because it is undetectable in the case of 
distant objects characterized by extremely broad interstellar Ly-a 
absorption lines that hide the heliospheric absorption. We know 
now that in the case of nearby objects with small interstellar col- 
umn densities, it can be detected in a number of directions. The 
heliospheric absorption will be very broad, thanks to the high tem- 
perature of the heliospheric H, and its centroid will also generally 
be shifted away from that of the interstellar absorption due to the 
deceleration, allowing its presence to be detected despite the fact 
that it remains blended with the interstellar absorption. The hydro- 
gen wall absorption was first detected by Linsky and Wood [1996] 
in Ly-a absorption spectra of the very nearby star a Cen taken by 
the Goddard High Resolution Spectrograph (GHRS) instrument on 
board the Hubble Space Telescope (HST). Since that time, it has 
been realized that the absorption can serve as a remote diagnostic 
of the heliospheric interface, and for stars in general, their "astro- 
spheric" interfaces. 
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models with the data, it was found that the kinetic models predict 
too much absorption. Models created assuming different values of 
the interstellar temperature and proton density fail to improve the 
agreement. Surprisingly, it was found that a model that uses a mul- 
tifluid treatment of the neutrals rather than the Boltzmann particle 
code is more consistent with the data [Wood et al, 2000b]. 

In this paper, we continue to study heliospheric absorption to- 
ward these stars. In our study, we use the Baranov-Malama model 
[Baranov and Malama, 1993, 1995, 1996] of the heliospheric in- 
terface, which is described in the next section. This model uses a 
Monte Carlo code with splitting of trajectories, which allows very 
precise computation of H atom distributions. Another advantage of 
the model is the possibility of separating of heliospheric H atoms 
into several populations, as discussed above. This model advantage 
allows us to consider separately two types of heliospheric absorp- 
tion, hydrogen wall absorption and heliosheath absorption. 

2. Heliospheric interface model 

To calculate heliospheric absorption we employ the Baranov- 
Malama model of the solar wind interaction with the two-component 
interstellar medium [Baranov and Malama, 1993, 1995, 1996; Iz- 
modenov, 2000]. This is an axisymmetric model, where the in- 
terstellar wind is assumed to have uniform parallel flow, and the 



Figure 1. The heliospheric interface is the region of the solar wind 
interaction with LIC. The heliopause is a contact discontinuity, 
which separates the plasma wind from interstellar plasmas. The 
termination shock decelerates the supersonic solar wind. The bow 
shock may also exist in the interstellar medium. The heliospheric 
interface can be divided into four regions with significantly different 
plasma properties: 1) supersonic solar wind; 2) subsonic solar wind 
in the region between the heliopause and termination shock (i.e., the 
heliosheath) ; 3) disturbed interstellar plasma region (or "pile-up" 
region) around the heliopause; 4) undisturbed interstellar medium. 



The a Cen line of sight lies 52° from the upwind direction of the 
interstellar flow through the heliosphere. An additional detection 
of heliospheric H I absorption only 12° from the upwind direction 
was provided by HST observations of 36 Oph [ Wood etai, 2000a] . 
For downwind lines of sight, the heliospheric absorption is also not 
negligible [Williams et al., 1997]. Analysis of Ly-a absorption to- 
ward Sirius shows that absorption of atoms of population 2 created 
in the heliosheath (see Figure 1) is needed in addition to interstel- 
lar and "hydrogen wall" absorption components in order to explain 
observations [Izmodenov et al, 1999]. 

In addition to heliospheric absorption, "astrospheric" absorption 
towards Sun-like stars was detected by Wood et al. [1996], Bring 
et al. [1997], and Wood and Linsky [1998]. These observations 
might be used to infer properties of these stars and their interstellar 
environments [Wood et al. 2001; Muller et al, 2001]. 

A theoretical model of the heliospheric interface should be em- 
ployed to interpret observations and put constraints on the local in- 
terstellar parameters and the heliospheric interface structure. This 
model should self-consistently take into account plasma and H atom 
components. Since the mean free path of H atoms is comparable 
with the size of the heliospheric interface, the H atom flow needs 
to be treated kinetically through the velocity distribution function, 
which is not Maxwellian [Baranov etal, 1998; Miiller et al, 2000; 
Izmodenov et al.,2QQV\. Recently, Woorf a/. [2000b] compared 
observations of Ly-a toward six nearby stars with model-predicted 
absorption. Both a Boltzmann mesh code [Lipatov et al, 1998] 
and a multi-fluid approach [Zank et al, 1996] were used to com- 
pute H atom distributions in the heliosphere. In comparing these 
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Figure 2. Density distributions of H atoms of population 3, atoms 
originated in the disturbed interstellar medium (region 3 in figure 
1), and population 2, atoms originated in the heliosheath (region 2 
in figure 2). Density distributions are shown as functions of helio- 
spheric orientation relative to different line of sights. Curves 1, 2, 3, 
4, 5, and 6 correspond to e = 12°, 52°, 73°, 112°, 139°, and 148°, 
respectively, where Q is the angle relative to the upwind direction 
of the interstellar flow. The figure shows the number densities for 
model 3 (see Table 1). 
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solar wind is assumed to be spherically symmetric at the Earth's 
orbit. Plasma and neutral components interact mainly by charge 
exchange. However, photoionization, solar gravity, and solar ra- 
diation pressure, which are especially important in the vicinity of 
the Sun, are also taken into account in the model. The process 
of electron impact ionization, which is important in the heliosheath 
[Baranov and Malama, 1996], is taken into account as well. Kinetic 
and hydrodynamic approaches were used for the neutral and plasma 
components, respectively. The kinetic equation for the neutrals is 
solved together with the Euler equations for a one-fluid plasma. 
The influence of the interstellar neutrals is taken into account in the 
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Figure 3. Number densities (upper row), velocities (middle row), 
and effective temperatures (lower row) of population 3 (left column) 
and population 2 (right column) toward 36 Oph. Dot-dash lines cor- 
respond to model 4, dot-dot-dot-dash lines correspond to model 5 
and long dashes correspond to model 6 (see Table 1). 

Table 1. Sets of model parameters'" 



Model 


nH,LlC 


"■p,LIC 


Notation in figures 


1 


0.10 


0.10 


solid 


2 


0.15 


0.05 


dotted 


3 


0.15 


0.10 


dashed 


4 


0.20 


0.05 


dot-dash 


5 


0.20 


0.10 


dot-dot-dot-dash 


6 


0.20 


0.20 


long dash 



" nH,LlC 2nd Up Lie are local interstellar atom and proton number 
densities, respectively. 



right-hand side of the Euler equations that contains source terms 
[Baranov and Malama, 1993; Malama, 1991], which are integrals 
of the H atom distribution function fHiVn) and can be calculated 
directly by a Monte Carlo method [Malama, 1991]. The set of 
kinetic and Euler equations is solved by an iterative procedure sug- 
gested by i?aranov e? a/. [1991]. Supersonic boundary conditions 
were used for the unperturbed interstellar plasma and for the solar 
wind plasma at Earth's orbit. The velocity distribution of interstellar 
atoms is assumed to be Maxwellian in the unperturbed LIC. 

Basic results of the model can be briefly summarized as follows 
[see also Izmodenov, 2000]. In the presence of interstellar H atoms 
the heliospheric interface is much closer to the Sun compared with 
the case of a fully ionized LIC. The termination shock becomes 
more spherical, and complicated shock structure in the tail disap- 
pears. The plasma flows are disturbed upstream of both the bow 
and termination shocks. This is an effect of H atoms that interact 
with protons by charge exchange. The solar wind is decelerated 
by 15 — 30% as it approaches the termination shock. The Voyager 
spacecraft observe the deceleration of the solar wind at large helio- 
centric distances [Richardson, 2001]. The number density of pickup 
ions, which are created by charge exchange, may reach 20 — 50% 
of the solar wind number density at the distance of the TS. The 
expected distance to the termination shock is 80 — 100 AU and de- 
pends on local interstellar parameters. The interstellar atoms are 
significantly disturbed in the interface. Since the velocities of new 
atoms created by charge exchange depend on properties of the local 
plasma, it is convenient to distinguish four populations of atoms 
depending on the place of origin, as described above. The velocity 
distribution function of H atoms can be represented as a sum of the 
distribution functions of these populations. Izmodenov et al. [200 1 ] 
studied the evolution of the four velocity distribution functions in 
the heliospheric interface. 

We computed the global heliospheric interface structure and dis- 
tribution of interstellar hydrogen in the interface using the Baranov- 
Malama model. Our calculations assume the following values for 
the fully ionized solar wind at 1 AU: rip, E — 6.5cm~'', V^,e — 450 
km s^^, and Me ~ 10, where np,E, V^,e and Me are the proton 
number density, solar wind speed, and Mach number, respectively. 
For the inflowing partially ionized interstellar gas, we assume a ve- 
locity of 25.6 km s^^ and a temperature of 7000 K. These values 
are consistent with in situ observations of interstellar helium [Witte 
et al., 1996]. Two other input parameters, interstellar proton and H 
atom number densities, are varied. The values assumed for these 
parameters for the various models are listed in Table 1. We vary 
Up, Lie in the range 0.1 — 0.2 cm^'^, while uh.lic is varied in the 
range 0.05 - 0.2 cm'^. 

In principle, each of the four H atom populations produces an 
absorption line. However, populations 1 and 4 can be neglected 
here, for the following reasons. Population 4 atoms are unperturbed 
interstellar atoms (having not charge exchanged at any time), and 
the average temperature and velocity of this population in the he- 
liosphere is very similar to those of the gas in the interstellar cloud 
beyond the heliosphere (there is a slight and negligible difference 
due to selection effects). Thus, the population 4 absorption cannot 
be distinguished from the interstellar absorption. Since the interstel- 
lar column-density is a fitted parameter, it will automatically include 
population 4. Moreover, this gas is not hot and the column density is 
small, resulting in an unnoticeable difference. Population 1 corre- 
sponds to atoms flowing away from the Sun at very large velocities 
(i.e., the speed of the neutralized supersonic solar wind protons). 
The produced absorption is thus strongly shifted and broadened due 
to solar wind inhomogeneities. Also, since column densities are 
very small, this absorption is very shallow and undetectable. 

Figure 2 shows calculated number densities of populations 2 and 
3, which are of interest in this paper, toward six selected line of sights 
from upwind to downwind. The figure shows the number densities 
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for model 3 (see Table 1). The hydrogen wall (upper panel) is most 
pronounced in the upwind direction. However, the wall decreases 
slowly with increasing of the angle 9 relative to the upwind direction 
{6 = 0° for upwind). MO ^ 112° (curve 4) the wall is still only 
slightly lower than it is in the upwind direction. There is no hydro- 
gen wall in downwind directions (curves 5 and 6). In upwind and 
crosswind directions, the number densities of population 2 are about 
two orders of magnitude less than number densities of population 3. 
However, the densities of populations 2 and 3 are different by only 
about one order of magnitude in the downwind directions (curves 5 
and 6). Despite their small densities, the atoms of population 2 can 
produce noticeable absorption due to their higher temperature. 

Figure 3 shows number densities, bulk velocities, and effective 
kinetic temperatures for models 4-6 toward 36 Oph. The highest 
hydrogen wall is for the model with the largest Up^Lic (model 6). 
Decrease of Up^Lic results in a decrease of the hydrogen wall. At 
the same time, the hydrogen wall becomes wider, so that the H col- 
umn density does not vary significantly from model 4 to model 6. 
All models show approximately the same velocities and tempera- 
tures of population 3 in the hydrogen wall region, which produces 
most of the heliospheric absorption. In contrast to population 3, the 
number densities, velocities, and effective kinetic temperatures of 
population 2 vary significantly with Up.Lic (right column plots in 
Figure 3). These effects are connected with: a) the location of the 
termination shock, which is closer to the Sun for higher Up^Lic, b) 
the width of the heliosheath region; and c) effects of interstellar H 
atoms filtration and their influence on the heliosheath plasma flow 
[for details, see Izmodenov, 2000; Izmodenov et al, 1999, 2001]. 

More results on the distribution of H atom and plasma 
parameters for the presented six models are provided at 
http://izmod.ipmnet.ru/ izmod/Papers/IWL/index.html. This page 
contains contour plots of plasma and H atom population number 
densities and velocities, and other supplemental information. 

3. Modeling of heliospheric absorption 

The absorption profile along a line of sight is 

7(A) = /o(A)exp(-r(A)), (1) 

where Io{X) is the assumed background Ly-a profile, and t(A) is 
the opacity profile. The opacity profile is 



Voigt profile: 



0(A). 



(2) 



Here, / is the oscillator absorption strength, A'^ is the column density, 
and (f>{\) = {XoC(f>{v)) / , where (f>{v) is the normalized velocity 
distribution along the line of sight. In the velocity space to which 
wavelengths can be linearly transferred, v — c(l — Aq/A), where 
Ao is the rest wavelength of Ly-a. In the case of a Maxwellian 
velocity distribution of the neutral gas, the absorption is determined 
by three parameters (A*', u, b), and the absorption profile is 



/(«) = Io{v)exp{' 



■KC^fNXo e 



" hi 



(3) 



Here, u and b are line-of-sight bulk and thermal velocities, respec- 
tively. Since the column density of H atoms in the heliospheric 
interface is less than 10^^ cm^'^, the extended Lorentzian wings of 
natural line broadening do not accumulate any appreciable opacity. 
Therefore, equation (3), which corresponds to pure Doppler broad- 
ening, is applicable. Note that for the interstellar absorption the 
column density is higher and one must take the Lorentzian wings 
into account and replace the Doppler profile in equation (3) with a 



I(v) = Io{v)exp{' 



mc^/TTb 



H{a, V — u)). 



(4) 



Here a — r/47r6, where F = 0.0755 is the transition rate in km/s 
units, and 



H{x,y) = 



exp{—t^)dt 



t)^+x^' 



(5) 



Using the distributions of interstellar H atoms computed as dis- 
cussed in the previous section, we compute the heliospheric absorp- 
tion toward the six selected directions. To study common properties 
and variations of the heliospheric absorption from upwind to down- 
wind, we assume in this section stellar profiles of /o(A) — 1. We 
compute absorption produced by H atoms of population 3 and pop- 
ulation 2 separately. Hereafter, these two absorption components 
will be called "hydrogen wall" and "heliosheath" absorption, re- 
spectively. 
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Figure 4. Simulated interstellar and heliospheric absorption of pop- 
ulations 2 and 3 toward 36 Oph (upper panel), 31 Com (middle 
panel), and Sirius (lower panel), which are located in the upwind 
(12°), crosswind (73°) and downwind (139°) directions, respec- 
tively. Interstellar absorption are shown by solid curves. Simulated 
hydrogen wall and heliosheath absorption are shown as dot-dash 
lines for model 4, dot-dot-dot-dash lines for model 5 and dash lines 
for model 6 (see Table 1). (The hydrogen wall absorption is the 
narrow, saturated absorption, while the heliosheath is the broader, 
unsaturated absorption.) 
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Figures 4 and 5 present both calculated hydrogen wall and he- 
liosheath absorption components for different models of the he- 
liospheric interface toward three stars, representing nearly upwind 
(36 Oph), crosswind (31 Com), and downwind (e Eri) directions. 
The absorption is shown in the heliospheric rest frame on a veloc- 
ity scale. The expected interstellar absorption toward these stars is 
also shown in the figures in order to compare them with the helio- 
spheric absorption. Details of the computation of the absorption are 
discussed in the next section. Figure 4 shows variations of helio- 
spheric absorption with increase of the interstellar proton number 
density, Up^Lic, from 0.05 cm^'^ (model 4) to 0.2 cm^"^ (model 6), 
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Figure 5. Same as figure 4, but for models 1 (solid lines), 3 (dashed 
lines), and 5 (dot-dot-dot-dashed lines). 

Table 2. Information on stars and interstellar absorption. 



Star Angle Distance Interstellar parameters toward the star" 





from 


(pc) 


u 


T 


N 


D/H, 




Upwind 




(kms-i) 


(K) 


X 10^^ 


X 10- 












(cm-2) 




36 Oph 


12° 


5.5 


-28.4 


6200 


7.80 


1.25 


a Cen 


52° 


1.34 


-18 


5800 


4.00 


1.5 


31 Com 


73° 


90.9 


-3.2 


8000 


8.3 


1.75 


/3Cas 


112° 


13.9 


9.1 


9500 


13.5 


1.65 


Sirius 


139° 


2.7 


19 


6000 


1.6 


1.65 








13 


6000 


1.6 


1.65 


e Eri 


148° 


3.3 


20.8 


7300 


7.6 


1.45 



" «,T,N are the interstellar velocity, temperature, column density and 
D/H ratio toward the star. 



with nH,Lic ~ 0.2 cm^ for all these models. Figure 5 shows 
variations of heliospheric absorption with an increase of nH,Lic 
from 0.10 cm"^ (model 1) through 0.15 cm"^ (model 3) to 0.2 
cm^^ (model 5), with Up.Lic ~ 0-1 cm~^ for all these models. 

Absorption of the hydrogen wall. All models produce nearly 
the same hydrogen wall absorption for upwind and crosswind direc- 
tions. Downwind, the absorption is different for different models. 
But due to small column densities of population 3 in downwind di- 
rections, the absorption cannot be detected in stellar spectra. In the 
upwind direction, the atoms of population 3 are decelerated relative 
to the original interstellar atoms. Therefore, the absorption by pop- 
ulation 3 is redshifted compared to the interstellar absorption. In 
nearly crosswind directions the projections of the bulk velocities to 
the direction of the Sun are close to zero for population 3 as well as 
for primary interstellar gas. Thus, the absorption of population 3 is 
hidden by saturated interstellar absorption and cannot be observed 
(middle panels in Figures 4 and 5). 

Absorption in the heliosheath. In contrast to the hydrogen wall 
absorption, the heliosheath absorption varies significantly with the 
interstellar proton and H atom number densities. A higher number 
density of population 2 for models with smaller interstellar proton 
number density Up^Lic (see Figure 3) leads to large absorption for 
this model compared to models with larger Up^Lic- Models with 
higher uh.lic also predict more heliosheath absorption (see Figure 
5) than models with small uh.lic- The difference between models 
increases for crosswind and downwind directions. Generally, for all 
models the heliosheath absorption is more pronounced in the cross- 
wind and downwind directions. Toward these directions both the 
size of the heliosheath and the number density of population 2 are 
larger compared with upwind. The heliosheath absorption is red- 
shifted in crosswind directions compared with the interstellar and 
hydrogen wall absorption components, and analysis of absorption 
toward crosswind stars might put particularly strong constraints on 
the heliosheath plasma structure. 

As seen from Figures 2 and 4, the parameters of heliospheric 
populations of H atoms have large variations in the heliospheric in- 
terface, as well as in each of its sub-regions. Moreover, the velocity 
distributions of the heliospheric populations are not Maxwellian [Iz- 
modenov et al. , 2000] . However, it is still interesting to see whether 
it is possible to model heliospheric absorption as if it were pro- 
duced by uniform Maxwellian gases. Using the computed distribu- 
tion functions, we take the zeroth, first, and second moments of the 
Maxwellian to compute A'^, u, and b. Then, we use these numbers 
to compute a Gaussian profile according to eq. (3). A comparison 
of the original profiles computed directly from the velocity distri- 
butions with the Gaussian profiles shows rather good agreement 
for both populations 2 and 3 (see Figure 6). Note that we com- 
pute Gaussian profiles and compare them with the exact profiles 
for populations 2 and 3 separately. If populations 2 and 3 were 
treated as one population, the Gaussian approximation would be 
much worse. Note that for some directions and models the differ- 
ence between model predictions and their Gaussian approximations 
might be larger, and these fits should be used for interpretation with 
great caution. 

It is important to note here that our model does not extend far 
enough downwind. Our computational grid extends only 700 AU in 
the downwind direction, which is not far enough to incorporate all 
heliospheric absorption [Izmodenov andAlexashov, in preparation]. 
Therefore, our calculations underestimate absorption in downwind 
directions (9 > 125°). 

4. The observations, interstellar absorption, and 
stellar spectra 

In this section, we use the heliospheric absorption models dis- 
cussed above to analyze absorption spectra toward six nearby stars. 
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For accurate study of the heliospheric absorption based on the anal- 
ysis of observed absorption spectra, one needs to know both the 
original stellar Ly-a line profile and the interstellar absorption. 

The Ly-Q lines of these six stars observed by HST are shown 
in Figure 7 (left column). All of the spectra except that of 36 
Oph A were taken using the Goddard High Resolution Spectro- 
graph (GHRS) instrument. The 36 Oph A data were obtained by 
the Space Telescope Imaging Spectrograph (STIS), which replaced 
GHRS in 1997. The spectra are plotted on a heliocentric velocity 
scale. All show broad, saturated H I absorption near line center and 
narrower deuterium (D I) absorption about 80 km s~^ blueward of 
the H I absorption. Figure 7 also shows the assumed intrinsic stellar 
Ly a lines and the best estimates for interstellar absorption. The 
interstellar parameters of these estimates are summarized in Table 
2. The models suggest that heliospheric absorption will always be 
redshifted from the the ISM absorption regardless of the direction 
of the line of sight, while astrospheric absorption will conversely 
be blueshifted [Wood et al, 2000b; Izmodenov et al. 1999]. This 
is primarily because of the deceleration and deflection of neutral H 
at the bow shock, and because for the heliosphere we are observing 
the decelerated and deflected neutral H from inside the heliosphere 
while for astrospheres we are observing it from the opposite per- 
spective outside the astrosphere. Thus, excess absorption on the red 
side of the Ly a absorption, when present, is interpreted as helio- 
spheric absorption, while excess absorption on the blue side, when 
present, is interpreted as astrospheric absorption. In this paper we 
focus on the red side of the Ly a absorption since we are interested 
in the heliospheric absorption. The stellar profiles and interstellar 
absorption estimates are based on previously published work, which 
is summarized below. 

a Cen. The a Cen data provided the first evidence for helio- 
spheric absorption. Both members of the ct Cen binary system were 
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from distributions of H atoms (solid lines) with absorption pro- 
file approximations assuming uniform Maxwellian gases (dashed 
lines). The absorption profiles are shown toward the 36 Oph direc- 
tion {6 = 12°). Upper and bottom panels are for populations 2 and 
3, respectively. 



observed. Linsky and Wood [1996] demonstrated that the results of 
their analysis were the same for both the a Cen A and a Cen B data. 
It is the a Cen B spectrum that we display in Figure 7, since it has a 
somewhat higher signal-to-noise ratio (S/N) than the a Cen A data. 
Using the D I line to define the central velocity and temperature of 
the interstellar material, Linsky and Wood [1996] found they could 
not fit the a Cen H I profile with only interstellar absorption. In 
Figure 7, we show interstellar absorption toward a Cen with the 
D/H value equal tol.5x 10^^, which is the generally accepted LIC 
D/H value. 

36 Oph. Wood et al. [2000a] analyzed the 36 Oph data. The 
analysis of the 36 Oph data proved to be very similar to that of a 
Cen, with the apparent existence of excess Ly-a absorption on both 
the red and blue sides of the line. Wood et al. [2000a] argue that 
only contributions of both heliospheric and astrospheric absorption 
can account for all of the excess absorption. 

Sirius The Sirius data were first presented by Benin et al. [ 
1995a, 1995b]. Unlike the other lines of sight presented here, there 
are two interstellar components instead of just one. Nevertheless, 
Berlin et al. [1995a] found that if they constrained the H I ab- 
sorption by assuming the LISM temperature and D/H ratio reported 
by Linsky et al. [1993] toward Capella (T = 7000 K and D/H = 
1.65 X 10^""), they could not account for excess H I absorption on 
both the blue and red sides of the line. Izmodenov et al. [1999] 
proposed that the red excess is due to heliospheric absorption and 
the blue excess is due to astrospheric absorption. 

31 Com, p Cas and e Eri. The 3 1 Com, /3 Cas and e Eri data 
were first analyzed by Dring et al. [1997]. The 31 Com and l3 Cas 
spectra can be fitted by single interstellar absorption components, 
but the e Eri spectrum is more complex. The H I absorption is 
significantly blueshifted relative to the D I absorption, indicating a 
substantial amount of excess H I absorption on the blue side of the 
line. Dring et al. [1997] interpreted this excess absorption to be 
astrospheric in origin. 

5. Comparison of theory and observations 

The middle and right columns of Figure 7 show the Ly a lines 
zoomed in on the red side of the HI absorption line, since that is 
where most of the heliospheric absorption is expected. The Ly a 
profile after interstellar absorption is shown in the plots as solid 
lines. The middle column also shows Ly a profiles after interstellar 
and hydrogen wall absorption, while the right column shows Ly a 
profiles after interstellar, hydrogen wall, and heliosheath absorption. 
The hydrogen wall and heliosheath components are shown for all 
six models discussed above (see Table 1). 

The absorption spectra toward upwind directions (36 Oph and 
a Cen) are fitted rather well by hydrogen wall absorption only. In 
comparing the theoretical absorption with data, the most important 
place for the model to agree with the data is at the base of the ab- 
sorption, where it is very difficult (if not impossible) to correct any 
discrepancies by altering the assumed stellar Ly a profile [Wood et 
al., 2000b]. Toward 36 Oph, model 2 produces slightly more ab- 
sorption and model 6 produces slightly less absorption at its base 
compared to the data. For other models the differences between 
model predictions and HST data at the base can be reduced by small 
corrections to the assumed stellar profile and interstellar absorption. 
Toward a Cen, all models show slightly more absorption at the base 
than the data, except for model 6. It is interesting to note that model 
6 fits the a Cen data remarkably well, while it does not fit as well 
the Ly a spectrum of 36 Oph. Note that these discrepancies with 
the data just discussed are small enough that they could in principle 
be corrected by reasonably small changes to the assumed stellar Ly 
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Figure 7. Left column: HST Ly a spectra of six stars: 36 Oph, a Cen, 3 1 Com, (5 Cas, Sirius, e Eri, respectively 
from top to bottom. Each plot shows the observed profile (thick solid line), assumed stellar line profile and interstel- 
lar absorption (thin solid lines). Middle column: Reproduction of left column, zoomed in on the red side of Ly-a 
absorption line. In addition, simulated Ly-a profiles after interstellar and hydrogen wall absorption are shown for 
models 1-6. Solid lines, when it is different from assumed stellar profile and interstellar absorption, correspond to 
model 1; dotted lines show results of model 2; dashes correspond to model 3; dot-dash curves correspond to model 4, 
dot-dot-dot-dash lines correspond to model 5, long dashes correspond to model 6. Right column: same as the middle 
column, but heliosheath absorption is added to the simulated Ly-a profiles. 
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too much absorption toward both 36 Oph and a Cen and produces 
the greatest difference with the data at the base of the absorption 
where it matters most. Consideration of the heliosheath absorption 
significantly increases the total heUospheric absorption predicted 
by model 4. The models with the worst agreement in the upwind 
directions (model 2 and model 4) correspond to the extreme case of 
small interstellar proton number density, Up^Lic- 

The observed absorption in the crosswind lines of sight to 31 
Com and /3 Cas shows no need for any heliospheric absorption, 
and our models do not generally predict significant absorption in 
the crosswind directions that can be clearly distinguished from in- 
terstellar absorption. Significant heliosheath absorption in model 4 
makes the model the worst fit to the data in the crosswind directions. 
The other models fit the data acceptably well. 

Results of the comparison for downwind lines of sight are more 
puzzling. It is clearly seen that the missing absorption toward Sirius 
can be easily explained by absorption in the heliosheath. Moreover, 
the heliosheath absorption predicted by different models is notice- 
ably different. Model 4, which has the worst fit to the data for all 
other lines of sight, fits the observed Sirius spectrum better than 
other models. 

In contrast to Sirius, for e Eri the interstellar absorption fits the 
observed spectrum very well and there is no need for heliospheric 
absorption. There is no hydrogen wall absorption in this direction 
(middle column of Figure 7), but all of our models predict too much 
heliosheath absorption. A similar problem was also found by Wood 
et al. [2000b] when they compared their kinetic models with the 
data. The discrepancy may be even more dramatic than our models 
suggest, because our computational grid extends only 700 AU in 
the downwind direction, which as stated above is not far enough 
to accoimt for all the absorption. A larger grid would presumably 
make the e Eri discrepancies worse for all models, and also change 
which model fits best for Sirius. 

6. Discussion 

As seen from the previous section, absorption spectra toward 
five of six stars can be explained by taking into account the hydro- 
gen wall absorption only. Only the Sirius line of sight requires a 
detectable amount of heliosheath absorption to fit the data. Unfor- 
tunately, the hydrogen wall absorption is not very sensitive to such 
interstellar parameters as the interstellar proton and H atom num- 
ber densities. The hydrogen wall absorption is most apparent and 
detectable in upwind directions. The small differences between the 
models and observations in these directions can be eliminated by 
small alterations of the stellar profiles (see Figure 7). The poorest 
agreement is for model 2. After including heliosheath absorption, 
model 4 seems to predict too much heliosheath absorption compared 
with the upwind data. Therefore, model 2 and model 4 produce the 
worst agreement with the data. This conclusion is also consistent 
with the analysis of absorption in crosswind directions (31 Com, (3 
Cas). 

Downwind lines of sight (Sirius and e Eri) can serve as good 
diagnostics of the heliosheath absorption. For the assumed stel- 
lar line profiles, all models predict too much absorption toward the 
most downwind line of sight, e Eri. However, since most of the 
discrepancy is away from the base of the absorption, we can try to 
correct this problem by modifying the shape of the assimied stellar 
Ly-a profile. Figure 8 demonstrates how the stellar profile of e Eri 
would have to be modified in order to improve agreement with the 
data. We do not display the central part of the profile, because the 
absorption is saturated, providing no information on the stellar Ly-a 
profile there. Because the models predict too much absorption, the 
fluxes of the assumed stellar profile must be increased to improve 
the fit. 

The revised profiles in Figure 8 are plausible in that they generally 
do not contain extremely steep slopes or fine structure. However, it 



is questionable if the tallest of these profiles (for models 4 and 5) are 
truly realistic. The revised stellar profiles in Figure 8 would suggest 
that the stellar Ly a profile of e Eri is significantly taller and nar- 
rower than the line profiles of similar stars like a Cen B and 36 Oph 
(see Figure 7), and the Sun. There is some subjectivity in deciding 
what is plausible and what is not, but we would conclude from Fig- 
ure 8 that models 4 and 5 probably require unrealistic modification 
to the stellar line profile to fit the data, but the other four models 
that require less extreme alterations of the assumed stellar profile 
may be all right. In truth, however, the situation is actually worse 
than this, because as mentioned above our models do not extend 
far enough downwind and therefore underestimate the amount of 
heliosheath absorption. Since a larger grid size would significantly 
worsen agreement with the data, it is uncertain whether any of the 
six models is truly consistent with the e Eri data. Furthermore, since 
the kinetic models presented by Wood et al. [2000b] had the same 
problem with the e Eri line of sight, it is possible that inaccuracies 
in the physical assumptions involved in current kinetic models may 
be resulting in overpredictions of absorption in downwind direc- 
tions, although additional theoretical work is required to test this 
interpretation. 

The models that have the least disagreement with the e Eri line 
of sight are models 1 and 6. The cormnon feature of these models is 
that both predict small number densities of H atoms of population 
2. In model 1 this is due to a low assumed interstellar H density, 
while in model 6 this is due to the small size of the heliosheath. 
Model 6 in fact predicts the narrowest heliosheath region of our 
six models. However, the interstellar proton and H atom number 
densities of this model seem to be unrealistically high, because the 
model predicts that the termination shock should be at 70 AU in 
the upwind direction. As of August 2001, Voyager 1 was at 82 AU 
and Voyager 2 was at ~65 AU, and neither has yet encountered the 
termination shock. Due to its low assumed interstellar H density, 
model 1 predicts the smallest amount of H atoms penetrating through 
the heliopause into the heliosphere of our six models. The number 
density of H atoms at the termination shock is 0.056 cm~^ for this 
model. This value is very close to the number density derived from 
the observed deceleration of the solar wind due to interaction with 
interstellar H atoms. Comparison of the solar wind speed measured 
on Ulysses and Voyager shows a decrease of about 40 km s~^, or 
10% in radial speed near 60 AU. Wang and Richardson [2001] have 
shown that this speed decrease implies an interstellar neutral density 
at the termination shock of 0.05 cm~*, which is close to our model 
1 value. 

Consideration of the Sirius line of sight complicates matters fur- 
ther. Model 2 and especially model 4 fit the observed profile very 
well, once again focusing our attention at the base of the absorption 
(see Figure 7), while the other models predict too little absorption. 
Models 2 and 4 assume a small interstellar ionization fraction in the 
vicinity of Sun, and recent studies of the LlC ionization state are in 
favor of such small ionization [Lallement, 1999]. However, model 
4 does not fit the other lines of sight well, and it must be noted once 
again that the models will underestimate the amount of heliosheath 
absorption in downwind directions due to the limited extent of the 
current model grid. 

The Ly-a absorption profile of the Sirius spectrum has been a 
controversial subject for the last few years and will probably remain 
so. First interpreted as due to absorption by a stellar wind from 
Sirius A [Bertin et al, 1995a], the extra absorption on the blue side 
of the line has later been explained as resulting from an astrosphere 
around Sirius and shown to be compatible with a very crude model 
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of such an astrosphere. However, the comparison between the Sir- 
ius B and Sirius A spectra later obtained by Hebrard et al. [1999] 
has shown that a large fraction (if not the totality) of the extra ab- 
sorption is specific to Sirius A and must arise closer to the star than 
predicted by an astrospheric model, favoring again a stellar wind 
as the source of the absorption. On the other hand, the extra ab- 
sorption on the red side, first interpreted as due to interstellar hot 
gas [Benin et al, 1995b], has later been compared satisfactorily 
with heliospheric absorption by the heliosheath [Izmodenov et al., 
1999], an interpretation which appears plausible based on Figure 7. 
Unfortunately, the comparison with Sirius B does not help here to 
disentangle the sources of the absorption, since the Sirius B spec- 
trum shows a broader and deeper absorption on the red wing, which 
encompasses the Sirius A line, and is identified by Hebrard et al. 
[1999] as due to the intrinsic photospheric Ly-a profile. Hebrard et 
al. [1999] also note that the excess absorption on the red side of the 
line can be removed if one relaxes the assumption that both inter- 
stellar components detected toward Sirius have D/H = 1.6x10^^. 
The heliosheath or interstellar gas remain possible sources for the 
excess absorption seen towards Sirius A on the red side of the line. 
If the absorption is due to the heliosheath, this is the only line of 
sight in which heliosheath absorption is detected. 

We now suggest one possible cause of our difficulties with down- 
wind lines of sight. The process of pickup ion assimilation into the 
solar wind is a very complicated phenomenon. Currently, it is be- 
lieved that pickup ions form a separate population of protons, which 
is co-moving with the solar wind proton plasma. The process of 
energy exchange between the populations is quite slow and it is ex- 
pected that the relaxation time is large compared with the time that 
is needed for a proton to reach the termination shock distance. How 
the non-equilibrium distribution of pickup ions is affected by the ter- 
mination shock is still an open question, although some scenarios 
have been developed [e.g., Fichtner, 2001]. Our model considers 
the solar protons and pickup ion protons as one fluid. This approach 
is based on fundamental questions of mass, momentum, and energy 
conservation. Therefore, it is quite natural to expect that the model 
predicts the positions of the shocks and the heliopause reasonably 
well. At the same time, measurements of speed and temperature 
of the solar wind at large heliocentric distances by Voyager are in 
favor of the two populations being co-moving, but thermally dif- 
ferent. This suggests a different scenario of the heliosheath plasma 
flow than is assumed in our model. 

This new scenario would result in different distributions of H 
atoms created in the heliosheath, which could possibly solve our 
difficulties in modeling the heliospheric absorption in downwind 
directions. A priori, we cannot exclude a scenario that provides 
the greatest absorption at heliocentric angles of about 130 — 150° 
from upwind. For such a scenario, perhaps all six lines of sight 
in Figure 7 could be explained simultaneously. Calculating ener- 
getic neutral atom (ENAs) fluxes, Gruntman et al. [2001] consider 
different plasma scenarios of heliosheath plasma flow. One of the 
models considered in the paper suggests a maximum of ENA fluxes 
in directions of 130 — 150° from upwind. This model assumes 
that the pickup proton population is carried through the shock with- 
out thermalization, i.e., preserving its velocity distribution filling the 
sphere in velocity space. Since ENAs originate in the heliosheath, it 
supports the idea to consider different scenarios of the heliosheath 
plasma in order to interpret all six observable spectra simultane- 
ously. It remains to be seen whether such models could both re- 
produce the absorption seen toward Sirius and solve our problem of 
overpredicting absorption toward e Eri. 

Finally, we have to note that an axisymmetric model of the helio- 
spheric interface was used to calculate the heliospheric absorption. 
The actual heliosphere is not axisymmetric due to both the latitudi- 
nal variation of the solar wind and possible effects of the interstellar 
magnetic field. Note, however, that the upwind direction of the ISM 



flow is only about 8° from the ecliptic, so variations of solar wind 
properties with latitude are unlikely to produce major differences 
from the axisymmetric approximation. In any case, asymmetry in 
the heliosheath or hydrogen wall could result in the heliospheric ab- 
sorption being different from axisymmetric. A 3D model is needed 
to assess the inaccuracies introduced by the axisymmetric assump- 
tion, but unfortunately there is as yet no adequate 3D model includ- 
ing the interstellar H atom population that can address this issue. 
Such a model is currently being developed by us. 

At the same time, our results can be interpreted as indirect evi- 
dence that there is no strong deviation of the heliospheric interface 
from being axisymmetric, since all four upwind/crosswind lines of 
sight discussed here can be well interpreted on the base of axisym- 
metric model. Taking into account the discussion above, one may 
conclude that our difficulties with the downwind lines of sight have 
other interpretations. However, this statement needs to be verified 
by 3D modeling of the heliospheric interface and by studying the 
sensitivity of heliospheric absorption to the 3D structure. 

7. Summary 

We have compared H I Ly-Q absorption profiles toward six nearby 
stars observed by HST with theoretical profiles computed using 
Baranov-Malama model of the heliospheric interface with six dif- 
ferent sets of model parameters. Our results are summarized as 
follows: 

1. It has been shown that the absorption produced by the hy- 
drogen wall does not depend significantly on local interstellar H 
atom and proton number densities for upwind and crosswind direc- 
tions. In downwind directions the hydrogen wall absorption is sen- 
sitive to interstellar densities, but this absorption component is most 
easily detected in upwind directions. In crosswind and downwind 
directions the hydrogen wall absorption is hidden in the saturated 
interstellar absorption and cannot be observed. 

2. The heliosheath absorption varies significantly with inter- 
stellar proton and H atom number densities. For all models, the 
heliosheath absorption is more pronounced in crosswind and down- 
wind directions. The heliosheath absorption is redshifted in cross- 
wind directions compared with the interstellar and hydrogen wall 
absorption components. 

3. Comparison of computations and data shows that all available 
absorption spectra, except that of Sirius, can be explained by tak- 
ing into account the hydrogen wall absorption only. Considering 
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heliosheath absorption, we find that all models have a tendency to 
overpredict heliosheath absorption in downwind directions. Toward 
upwind and crosswind stars the small differences between model 
predictions and the data can be corrected by small alterations of the 
assumed stellar Ly-a profile. However, the downwind e Eri Une 
of sight is a problem, as the models predict too much heliosheath 
absorption in that direction, and for many, if not most, of the models 
the discrepancy with the data is too great to resolve by reasonable 
alterations of the stellar profile. 

4. It is puzzling that model 4 provides the best fit to the absorp- 
tion profile toward Sirius but the worst fit to the other lines of sight. 
This may be due to our models underestimating heliosheath absorp- 
tion in downwind directions due to limited grid size, or perhaps the 
detected excess absorption toward Sirius is not really heliospheric 
in origin, as suggested by previous authors. It is also possible that 
the difficulties the models have with the downwind lines of sight 
towards Sirius and e Eri might be resolved by modifications to the 
models, perhaps by taking into accoimt the multi-component nature 
of the heliosheath plasma flow. 
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